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Abstract 



We consider incorporation of particle detachment in Smoluchowski 
model of colloidal growth. Two approaches are considered, utilizing phe- 
nomenological rate equation and exact large-time results. Our main con- 
clusion is that the value of the large-time diffusing particle concentration 
at the aggregate surface is the only parameter needed to describe the 
added effect of detachment. Explicit expression is given for the particle 
intake rate. 
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1. Introduction 



In this work we address aspects of the theory of growth of aggregates 
by capture of small diffusing particles (singlets), such as atoms, molecules, 
etc., with focus on processes relevant in colloidal growth. The celebrated 
work of Smoluchowski [1], reviewed, e.g., in [2], considered the process 
of diffusional capture of particles as follows. The growing aggregate is 
modeled "adiabatically" as a sphere of radius R. The concentration of 
particles is initially c p (p stands for particles). Next, the radially symmet- 
ric concentration of particles at time t, c(r, t), is calculated by solving the 
diffusion equation with the absorbing boundary condition at r = R, i.e., 
c(R, t > 0) = 0. The resulting particle intake 3> (number per unit time) 
into the sphere is usually taken as the large-time (steady-state) limit of 
the following expression obtained from this solution: 



where D is the singlet diffusion constant. The value of $(oo) is then used 
to estimate the growth rate of the spherical aggregate. 

In colloidal growth, a popular model has been that of "burst nucle- 
ation" by LaMer and Dinegar [3-4]. In this model, one adopts concepts 
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from the classical nucleation theory to colloidal aggregation. Specifically, 
for sufficiently small aggregate radii R, it is assumed that the aggre- 
gate sizes are distributed according to the equilibrium condition with the 
Boltzmann factor determined by the aggregate free energy G(R). This 
equilibrium distribution is usually assumed to apply for R < R c , where 
the critical radius R c corresponds to the maximum of the free energy. 

Larger aggregates, those with R > R c , are assumed to grow irre- 
versibly by the diffusional mechanism (1.1) and possibly by cluster-cluster 
aggregation, depending on the precise nature of forces between aggre- 
gates. The whole picture is "adiabatic" in that many of the quantities 
involved, specifically, G, R c , c p , evolve in time, presumably on the time 
scales much larger than the times needed to establish the asymptotic 
large-time diffusional flux; see (1.1). Thus, for the diffusion part of the 
description, the steady state results nearly always suffice. 

For clusters with R < R c , there must be dynamical mechanisms that 
lead to equilibration, if the "nucleation" droplet-model concept is indeed 
applicable to colloids. In addition to the mechanisms that correspond to 
internal rearrangement and restructuring, there must be also detachment 
of particles. Indeed, the aggregates cannot only grow, by capturing sin- 
glets, while at the same time achieve equilibrium distribution. Since the 
critical cluster size R c can be actually quite large, it follows "by continu- 
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ity" that detachment will also occur for clusters with R > R c which are 
not in equilibrium. 

The main goal of the present work is to incorporate detachment pro- 
cess in Smoluchowski model. Thus, we focus on the diffusion problem. 
We assume that the aggregates are large, R > R c , and that their radii 
vary slowly with time so that the time scales of singlet diffusional pro- 
cesses are fast in comparison. Furthermore, by focusing on the diffusion 
problem we avoid discussion of cluster-cluster aggregation and of detach- 
ment of multiplets, as well as of internal restructuring processes. The 
latter are present not only for small clusters but also for large clusters 
because without internal restructuring, diffusion-limited colloidal growth 
would result in fractals, e.g. [5]. Under properly chosen experimental 
conditions, e.g. [6] and references therein, the aggregates are as dense 
as the bulk material and have definite shapes, frequently spherical as 
assumed here, and semicrystalline structure. 

Our second goal is to emphasize a difference between the colloidal 
growth and classical nucleation. Specifically, we note that the intake ex- 
pression (1.1) for t = oo is proportional to R rather than to the surface 
area of the collector aggregate AttR 2 . In ordinary metastable-to-stable 
two-phase nucleation, when the transport mechanisms are not diffusional, 
for large clusters one can assume that each unit of the surface area cap- 



-4- 



tures and emits equal amounts of atoms or molecules. The net intake rate 
is then proportional to the surface area [7,8]. We would like to explore 
failure of this "surface extensivity" argument for colloids and include de- 
tachment in this discussion. A related matter is the mere existence of 
the "surface layer" since intuitively it seems to be connected with the 
"per-unit-area" adsorption and detachment. 

Finally, we would like to have few-parameter approaches that retain 
the simplicity and exact solvability of Smoluchowski model. We note that 
numerous modifications of Smoluchowski model have been proposed in 
the literature; see [2,9] for reviews. The most studied case in the field 
of colloids has been deposition of particles in a low- or intermediate- 
density monolayer on a solid substrate, without incorporation of the de- 
posited particle into the material of the collector by internal restructur- 
ing. Detachment from such monolayers is approximately proportional 
to the amount of matter deposited. It is modeled by phenomenological 
rate equations for the surface coverage coupled with expressions for the 
diffusional or convective-diffusional flux [10,11]. A detailed mathemati- 
cal study of one such boundary condition was published by Agmon [12]. 
Experimental results have been reviewed in [11,13]. 

Another application of Smoluchowski approach has been in the field 
of chemical reaction rates in the diffusion-limited regime, reviewed in 
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[9]. Here R is no longer large. Instead of detachment one has partial 
reaction, i.e., some particles may not react on encountering a collector. 
This is mathematically described by the radiation boundary condition 
[14]. 

The outline of the rest of this article is as follows. In Section 2, 
we consider a phenomenological model for the boundary condition at the 
aggregate surface that allows for detachment in Smoluchowski model. In 
Section 3, we derive a generalized radiation boundary condition which 
is exact in the large-time limit. The diffusion problem with the latter 
boundary condition is solved in Section 4. Finally, in Section 5 we com- 
pare the two approaches and conclude that for all practical purposes one 
can actually use a simpler boundary condition which only involves one 
new parameter: the concentration of the diffusing particles right at the 
aggregate surface. 
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2. Phenomenological Rate Equation Model 

In the original Smoluchowski model [1], the radially symmetric ver- 
sion of the diffusion equation, 

d -l = r ->°- (r*D^) 
dt dr \ dr J 

is solved for r > R, with the initial condition 

c( r ,t = 0) = c p (2.2) 

and boundary condition 

c(r = R,t>0) = 0. (2.3) 

The discontinuity in the value of the concentration at r = R (at t = 
as compared to that at t > 0), shared by many modifications of the 
problem, leads to singular behavior at t = 0, see (1.1), that has been of 
some concern in the literature on chemical reaction applications [2,9]. 

One can consider various modifications of the relations (2.1)-(2.3) 
to improve the model and incorporate effects other than the irreversible 
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(2.1) 



adsorption at r = R expressed by (2.3). In this section we consider a phe- 
nomenological modification of the boundary condition (2.3) to allow for 
detachment. We use Langmuir-type mean-field rate equation approach 
popular in applications of this sort in colloid literature, e.g. [15]. Thus, 
we propose the relation 



to replace (2.3). Here we assume that diffusing particles that reach the 
surface are incorporated in the aggregate structure at the rate Kc pro- 
portional to their concentration at R. The second term corresponds to 
detachment, and we assume that this rate only depends on the internal 
processes so there is no dependence on the external diffuser concentration. 

In order to solve (2.1) with (2.2), (2.4), we define Laplace Transform 
of the concentration, 




at r = R 



(2.4) 




(2.5) 



which satisfies the equation 




(2.6) 



8 



Here we already used the initial condition. The general solution of this 
equation is given by 



„3 + M e -M + M e M, (2.7) 
s r r 

where we set B = to have a bounded solution as r — > oo. The boundary 
condition at R becomes 



sn — c p = —Kn + ks 1 , at r = R . (2.8) 
The final result for the concentration is 



n = * - ( *) ^± e -V^(r-R) , (2.9) 
s \r J s(s + K) v ' 

The full time-dependence of the concentration can be expressed in terms 
of error functions, etc., by using the tables of Inverse Laplace Transform. 
However, the result is quite cumbersome and unilluminating. We present 
a simpler limiting expression shortly. 

For Laplace Transform of the intake rate 3>(t), denoted (f)(s), we get 



= AirDR * £ (l + R^/TjD) , (2.10) 

where we used 
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= AttR 2 D 




r=R 



(2.11) 



The time-dependence of the intake rate is then 



$(t) =AttDR{c p K -k) 



1 - e 
K 



-Kt 



+ 



VttDK 



2R 



L 



dx e 



y 

[2.12) 



We note that this rate expression is actually regular for small times. 
However, for colloid applications we are interested in large times. For 
t ^> 1/K, the limiting form of (2.12) is quite similar to Smoluchowski 
result (1.1). Keeping only the leading time-dependent term, we have 



The only change is the reduction of the rate due to detachment, propor- 
tional to the ratio k/K (the intake rate can actually become negative if 
the detachment is fast enough). 

We now turn to the concentration and present its leading time de- 
pendence in the limit t^> 1/K, 




(2.13) 




(2.14) 
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As t — > oo, the concentration profile develops a long-range tail ~ R/r that 
has no length scale of order 1 associated with it. In fact, the perturbation 
of the diffuser particle concentration extends for distances comparable to 
R. The vanishing of the term ~ R 2 in (2.13), so that only the "nonex- 
tensive" contribution ~ R survives as t — > oo, is associated with this 
property. Obviously, it is dependent on the transport mechanism being 
diffusional. Note that the terms retained in (2.13), (2.14) vary on the 
time scales set by R 2 /D. Since we are interested in large R values, this 
quantity is much larger than the "transient" time scale 1/K present in 
the exact expressions. 
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3. Model with Boundary Layer 

In this section we address the question of whether the boundary 
condition at r = R can be derived from a diffusional model. We will 
keep the assumption that the diffusion is free for r > R in the sense 
that there is no external potential. However, let us assume that in a 
certain boundary region from R — t to R, where I « fi, the particles 
move in the potential energy U (r) owning to their interaction with the 
aggregate. As r — > R — £ from above, the diffusional picture breaks down 
as the particles are subject to the internal restructuring mechanisms in 
the aggregate. We will assume that this dense aggregate structure can 
be represented by the boundary condition 

c(R-l,t) = c . (3.1) 

In fact, the diffusion constant D will be affected by the increased particle 
concentration, proximity to the semicrystalline aggregate, and external 
potential which itself can be partly functionally c-dependent. We ignore 
all these complications and keep D constant. 

Thus, our model is presented in Figure 1. We assume that Co ^> c p 
is a fixed quantity that only depends on the aggregate material. We also 



put 



V(r) = U{r)/k B T 



(3.2) 



and assume that V(r) is a fixed function of the difference r — R, vanishing 
for r > R. This potential energy (in units of ksT) may represent a 
truncated version of the full interaction energy at r. For "free" diffusional 
aggregation, we can assume that the interactions for r > R are weak. 
However, in the boundary layer of size £, large negative energy is expected, 
as sketched in Figure 1. Thus, we denote 



with Vo ^> 1. Furthermore, we expect that |V(r)| ^ 1 in most of the 
boundary layer. 

The diffusion equation is replaced by one with the potential energy 
term: 



We now offer exact considerations which are valid only in the steady-state 
t — > oo limit. Then the diffusion equation can be integrated as follows: 



V(R - £) 



Vo 



(3.3) 




(3.4) 




4yrr 2 



(3.5) 



Here c(r) denotes c(r, oo), and 3> = $(oo) is the (constant) intake rate 
into a sphere of radius r at time t = oo. 

Now the potential V(r) is continuous but its derivative is likely dis- 
continuous; see Figure 1. It then follows from (3.5) that the concentra- 
tion c(r) is continuous but the derivative dc/dr is not. Consider the limit 
r — > R from above, denoted by + here. We have the following expression 
involving the derivative: 



Next, we solve the steady-state diffusion equation for r < R for the 
concentration c(r) by direct integration and use of the boundary values 
(3.1) and (3.3). We also use V(R) = to get, for the the concentration 
at R, 






(3.7) 



Finally, we utilize (3.6) to obtain the following relation: 
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c(R) = a 




r^R+ 



+ 



(3.8) 



where 



a = R 




R e V(r) 



(3.9) 



(3 = c e 



-v 



(3.10) 



Since relation (3.8) applies with the derivative at R+, it can be 
used as the boundary condition for the free diffusion at r > R. We can 
thus lump all the effects of the boundary layer in two parameters a and 
(3. Such boundary condition resembles the radiation boundary condition 
used in the chemical reaction applications [9,14], which assumes that the 
concentration is proportional to its derivative. However, there are two 
important differences. Firstly, there is the additional constant term (3. 
Secondly, the coefficient a in our case is only weakly dependent on R. 
The integral in (3.9) is obviously of order 



Thus, the main result of this section is that the problem with the bound- 
ary layer can be replaced, at least for large times, by an equivalent prob- 



er ~ £e~ Vo . 



(3.11) 



lem of free diffusion at r > it*, as long as we use our generalized radiation 
boundary condition at R, 



Oc 

c = a — +/? , &t r = R. (3.12) 
or 
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4. Generalized Radiation Boundary Condition 



Before comparing, in the next section, the two choices of the bound- 
ary condition: phenomenological (2.4) and radiation (3.12), let us work 
out, in this section, what would be the time dependence of the free diffu- 
sion problem at r > R if the boundary condition (3.12) were used for all 
times. We go through steps similar to Section 2 and only emphasize the 
final expressions. In Laplace Transform version, the boundary condition 
becomes 

n = a — + (3s~ 1 1 at r = R . (4.1) 
or 

For the concentration and intake rate, we get 



s \ r 



s(l + aR- 1 + oty/s/D^j 



o = 4ttDR — v (l + RyflfD) . (4.3) 

s (l + aR- 1 + Oiy/s/D) V J 



These results can be inverted to get the time dependence since they 
only involve tabulated functions. For instance, 



= ^7rDR(c p -f3) 



a + R 



R 



l + ^ e («+^) 2 «- 2 ^- 2 ^Erfc 



a 



a + R 
aR 



y/Dt 




This expression is well-behaved at t = 0. As in Section 2, though, it is 
more useful to analyze limiting results for large R and large times. We 
note that (3 is the product of a large quantity (co ^> c p ) and a small 
Boltzmann factor (e~ v ° <C 1); see (3.10). However, a is the product 
of two small quantities: I « and the Boltzmann factor; see (3.11). 
Therefore, we can always assume that R ^> a. Given this assumption, 
one then identifies the transient time scale via the condition 



t > a 2 /D . 



(4.5) 



The asymptotic expression for the intake rate becomes 




(4.6) 



which is the same as (2.13) provided we identify 



k 
K 



(4.7) 



Similarly, with (4.7), the concentration is given by (2.14). 



5. Constant Boundary Condition. Discussion 

In view of the above discussion it is tempting to combine the phe- 
nomenological model which was introduced for all t- values with the exact 
expression for the boundary condition at t = oo. We could write some- 
thing like 



Here (3 is given by (3.10), a by (3.9) and (3.11), and it is also tempting 
to assume that there is only one transient time scale so that 



However, as long as we are interested only in the large R and large 
time behavior, for all practical purposes we can use instead the boundary 
condition 




at r = R . 



(5.1) 



K (xD/a 2 . 



(5.2) 



c = /3 



at r = R . 



(5.3) 



Indeed, the value of the concentration at R for large times, which is 
exactly (3 = k/K, independent of R, in both models, see (2.9) and (4.2), is 
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the only parameter needed to calculate the modification of the asymptotic 
intake rate due to detachment. With this boundary condition for all 
times, the results (2.13) or (4.6), and (2.14), become exact. The price 
paid is that the singularity at t = is back. However, only large t results 
are of interest in colloidal growth applications. 

In summary, a combination of phenomenological and exact consid- 
erations yields expressions for the diffusional intake rate of particles in 
Smoluchowski model with added detachment. The intake rate is reduced 
(it can become negative) and its value can be calculated from the large- 
time concentration at the collector surface {(3) and the value far away from 
the surface (c p ). The only added parameter is thus (3. At large times, the 
concentration profile develops a slowly varying contribution ~ R/r and 
the perturbation of the concentration due to the collector surface extends 
on length scales ~ R, far beyond the dimension of the "boundary layer" 
f « i?. This added concentration variation is not present in classical 
nucleation models with nondiffusional transport. The main result is the 
nonextensive (~ R instead of ~ R 2 ) behavior of the intake rate in the 
t — > oo limit in the diffusional model. 
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Figure Caption 

Figure 1: Schematic representation of the truncated dimensionless 
potential energy V(r), defined in Section 3. 
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